Inverse flux quantum periodicity of magnetoresistance oscillations 
in two-dimensional short-period surface superlattices 

X. F. Wang,^ P. Vasilopoulos*, and F. M. Peeters* 



''"^Concordia University Department of Physics 
1455 de Maisonneuve Quest, 
Montreal, Quebec, Canada, H3G 1M8 

*Departement Natuurkunde Universiteit Antwerpen (UIA), 
Universiteitsplein 1, B-2610, Belgium 

(Dated: February 2, 2008) 

Abstract 

Transport properties of the two-dimensional electron gas (2DEG) are considered in the pres- 
ence of a perpendicular magnetic field B and of a weak two-dimensional (2D) periodic potential 
modulation in the 2DEG plane. The symmetry of the latter is rectangular or hexagonal. The well- 
known solution of the corresponding tight-binding equation shows that each Landau level splits 
into several subbands when a rational number of flux quanta h/e pierces the unit cell and that 
the corresponding gaps are exponentially small. Assuming the latter are closed due to disorder 
gives analytical wave functions and simplifies considerably the evaluation of the magnetoresistivity 
tensor p^j^^. The relative phase of the oscillations in p^x and pyy depends on the modulation periods 
involved. For a 2D modulation with a short period < 100 nm, in addition to the Weiss oscillations 
the collisional contribution to the conductivity and consequently the tensor p^jj show prominent 
peaks when one flux quantum h/e passes through an integral number of unit cells in good agreement 
with recent experiments. For periods 300 — 400 nm long used in early experiments, these peaks 
occur at fields 10 — 25 times smaller than those of the Weiss oscillations and are not resolved. 

PACS numbers: 73.20.At; 73.20.Dx; 73.61.-r 
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I. INTRODUCTION 



In the last decade the magnetotransport of the 2DEG, subjected to periodic potential 
modulations, has attracted considerable experimental [l| and theoretical attention. 
For one-dimensional (ID) modulations novel oscillations of the magnetoresistivity tensor p^i, 
have been observed, at low magnetic fields S, distinctly different in period and temperature 
dependence from the usual Shubnikov-de Haas (SdH) ones observed at higher B. These novel 
oscillations reflect the commensurability between two length scales: the cyclotron diameter 
at the Fermi level 2Rc = l^fTjmJ^ ^ where ne is the electron density, i the magnetic length, 
and a the period of the potential modulation. The situation is similar but less clearcut for 
2D modulations from both a theoretical and an experimental [lU^l^oint of view. To 




date most of the experimental results pertinent to 2D modulations [7|]-[9[ with square or 
hexagonal symmetry have indicated strongly that the predicted jl^ fine structure of the 
Landau levels is not resolved. Magnetotransport theories pertinent to this case are rather 
limited in contrast with those for ID modulations. 

Recent observations !9l call for additional theoretical work since they could not be fully 

n 

explained by earlier semiclassical theories jj]. In this paper we develop, along the lines of 
Ref. 3, the relevant quantum mechanical magnetotransport theory of the 2DEG for precisely 
the case that the fine structure of the Landau levels is not resolved. Our goal is to explain 
recent experimental results [14] on 2D, short-period (a ~ lOOOA) surface superlattices with 
mobility [i in the intermediate range, i.e., /i ~ 100 m^/Vs. The symmetry of the 2D 
modulation is taken to be rectangular or hexagonal. A brief semiclassical account, pertinent 
to the former symmetry, was reported in Ref. 4 b). New magnetoresistance oscillations are 
found to occur vjhen one flux quantum h/e passes through an integral number of unit cells as 
was recently observed experimentally IJ]. Here we show that these oscillations result from 
the interplay between band conduction and collisional conduction. A new contribution to 
the latter opens up as hopping between cyclotron orbits which are separated by an integral 
multiple of the modulation period and have the same position relative to the modulation 
lattice. This contribution is appreciable only in short-period superlattices and accordingly 
could not be resolved in early experiments on long-period superlattices. 

In the next section we derive the one-electron eigenfunctions and eigenvalues for rectan- 
gular and hexagonal modulations; we also present the density of states. The analytical and 
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numerical results for the corresponding conductivity or resistivity components are presented 
in Sec. III. Numerical results are given in Sec. IV and concluding remarks in Sec. V. 



II. EIGENVECTORS, EIGENVALUES, AND DENSITY OF STATES 

We consider a 2DEG, in the (x, y) plane, in the presence of a perpendicular magnetic field 
B = Bz and of a 2D periodic potential modulation U{x,y). The electrons are considered as 
free particles with an effective mass m*. In the absence of the modulation the normalized 
one-electron eigenfunction, in the gauge A = {0,Bx,0), is given by e*^"^0„(a; + x^j ^fl7y 
where 0n(a^ + a^o) is the well-known harmonic oscillator function, centered at — Xq = —(?ky, 
and Ly is the sample's width. 

In the presence of a sinusoidal ID modulation one can use perturbation theory to 
evaluate the energy spectrum and eigenf unctions. Alternatively, one can use a tight-binding 
scheme, along the lines of Ref. and look for solutions of the one-electron Hamiltonian 
= (p + eA)^/2m* + \4 cos{Kxx) that are linear combinations of the unperturbed (V^ = 0) 
ones: | ipnky >= J2p^p I + pG >, where G is a suitable wave vector introduced here 
for convenience and \n,ky + pG > is the unperturbed state {K^ = 2Tr/a^, and is the 
modulation period along x). As in Ref. we take G = Ky = 271 /ay with ay the modulation 
period along y. The summation over p has to be extended to all integer p values such that 
—Lx/21'^ < ky+pKy < Lx/2i'^, where is the length. For p = we have the limits for ky as 
—ax/2i'^ < ky < ax/2i'^. Then the tight-binding equation < n, ky+pKy \ — E\ipnky >= 0, 
in which mixing of different Landau levels n is neglected, gives acceptable solutions for the 
coefficients Ap as Ap = Aq exp(i^p). The new states are labeled with the additional quantum 
number ^ (0 < ^ < 27r): | ipnky^ >= Xlp ^^P(^^P)I^5 ^2/ + pKy >. The orthonormality 
condition < ipnkyS, \ '^nky^' >= S^^' gives C, = 2TTh'i'^Ky/Lx — > k^i'^Ky, u being an integer, and 
Aq = i(Ky/ LxY^"^ by normalization. The energy spectrum obtained in this way is the same 
as that obtained by perturbation theory 

We will now use this information to obtain the corresponding eigenvectors and eigenvalues 
for a 2D modulation potential. 
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A. Rectangular symmetry 



We assume the following one-electron Hamiltonian 

H° = {p + eA)/2m* + 14 cos{K^x) + Vy cos{Kyy) (1) 

where = 27r/a^(/i = x,y)] and ay are the periods along the x and y directions, 
respectively. 

In the gauge chosen, A = (0, -Bx, 0), the second term of Eq. (1) is not diagonal in ky 
and therefore \ n,ky > is not a convenient basis set. But as in the ID case we can look 
for solutions of Eq. (1) in the form | fn,ky >= J2p^p I n,ky + pKy > as described above. 
This choice of eigenf unctions is also suggested by the fact that Vy cos Kyy connects the 
unperturbed state | nky > with only the states | n,ky± Ky >. In this case the equation 
< n, ky + pKy I — E I ipn,ky >= 0, in which mixing of different Landau levels is neglected, 
takes the form 

V^F^M cos(27rpa + K^Xo)Ap + ^VyFn{uy){Ap+i + Ap^i) = {E ~ En)Ap, (2) 

where a = 2'n'i'^ /a^ay, E^ = {n + l/2)hujc is the "unperturbed" eigenvalue and uj^ = \e\B/m* 
the cyclotron frequency. Further, = exp(— m^/2)L„(u^), is the Laguerre 

polynomial, and = 

The solution of Eq. (2) gives the eigenvalues E and the eigenvectors Ap. We see imme- 
diately that for a integer the equation admits the exponential solutions Ap = Aqc*^^, with 
Aq and ^ given above. This is also the case for those values of a for which Fn{ux) vanishes 
since = 27r^£^/a^ = iT{ay/a^)a. In the former case we have 

Enki = En + VxFn{ux) cos{Kr,xo) + VyF^iuy) COS i (3) 

and in the latter 

Enki = En + VyFn{Uy) COS^, (4) 

where ^ = {2'kv / Lx)l'^Ky = i'^Kyk^. In both cases the unperturbed Landau levels broaden 
into bands (with a bandwidth equal to 2(14|F„(m2;)| -|- Vy\Fn{uy)\) and 2Vy\Fn{uy)\, respec- 
tively), that oscillates with magnetic field B and (large) index n, cf. Refs. 2 and 3. The 
energy spectrum given by Eq. (3), plotted in Fig. [Hfor n = 0, a = 1, = 2Vy = 1 meV, 
is a periodic function of kx and of ky since ^ = i'^Kykx and xq = i'^ky. Notice that the 
arguments of the cosines in Eq. (3) can be shifted by 2na, a integer. 
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FIG. 1: Energy spectrum as a function of the wave vectors kx and ky for n = 0, a = 1, Vx = 2Vy = 1 
meV. The modulation wave vectors Kx and Ky are shown by the thick arrows. 

One important consequence of this nonzero bandwidth is that the mean velocities Vx and 
Vy, which vanish in the absence of modulation, are now finite: from = dEnk^/hdkfj_, fi = 
X, y, we obtain 

Vx = -{eKyVy/h)Fn{Uy) SiU ^ , (S) 

and 

Vy = -{fKxVx/h)Fniux) sm{KxXo). (6) 

Eqs. (5) and (6) lead to a finite diffusion or band conductivity which is absent when the 
modulation is not present, cf. Sec. II. 

Equation (2) is the same as Harper's equation but the coefficients V^-F„(m^) depend on 



the magnetic fielc 
by Hofstadter 



nown 



ll|b, 



. For B values other than those pertaining to Eq. (3) it has been si 
ll|a for the case of constant coefficients and by Claro and Wannier 
for the case that the latter depend on B (hexagonal modulation), that the energy spectrum 
resulting from the numerical solution of Eq. (2), i.e., E, shows, when E is measured in units 
of Vfj,Fn{u^), a nontrivial structure: for a = i/j, i,j being integers, each Landau level is 
split into j subbands and Eq. (2) is periodic with period j. Here, in view of the reported 



experiments 



which did not indicate that this fine structure of the energy spectrum was 
resolved, we will assume that this is indeed the case, i.e., that in samples of not exceptionally 
high mobilities, such as those of Ref. 5, the small gaps mentioned above are closed due to 
disorder and justify the assumption below, see subsection C. That is, we assume that the 
Landau levels are bands. Now from the numerical solution of Eq. (2) we know that the 
bandwidth corresponding to a = i/j cannot exceed the value obtained from Eq. (3) or Eq. 
(4). Therefore, for computational convenience, we will assume that the energy spectrum is 
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a. 







FIG. 2: Energy spectrum, obtained from the exact solution of Eq. ((2}, in (a) and (c), as a function 
of the wave vectors kx and ky for n = 0. The corresponding approximate spectrum, given by Eq. 
(3), is shown in (b) and (d). In (a) and (b) we have a = 1/2 and a = 2/3, in (c) and (d) a = 2/3. 
T,.e pe.„.. a,, a,. «00A pe«ai„ .o «,e exp.i.en. of Ref. Q. 

given by Eq. (3) and the eigenfunctions by | ipnky^ >= J2p ^^^^ \ n,ky + pKy >. 

In Fig. El we compare the energy spectrum obtained by exactly solving Eq. with the 
one given by Eq. (3). We do so because in the conductivity calculations we will use Eq. (3), 
as an approximation that will be justified, for all magnetic fields or values of a. The exact 
spectrum ior a = i/j is composed of j minibands. It's dependence on k^, which does not 
appear in Eq. (2), is obtained by introducing appropriate new basis states, in the manner 
of Ref. Q|, with \kx\ < n/iax and \ky\ < n/ay restricted in the magnetic Brillouin zone. 
As can be seen, the two spectra are quite different from each other. The corresponding 
difference in the density of states is much weaker if a small broadening is included, see 
subsection C below. For a integer, however, the exact spectrum and that given by Eq. (3) 
coincide; the result is shown in Fig. 1. 

There are two alternative, approximate ways to obtain Eq. (3). First, we take ~ 
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and use the corresponding ID tight-binding states \ n,ky,^ > to obtain the energy spectrum 
given by En + VxFn{ux) cos^K^xo). We then use first-order perturbation theory, involving the 
states I nky$, >, to evaluate the energy correction to this spectrum due to the term Vy cos Kyy 
for Vy <^ hjJc + Vx] the result is identical with that given by Eq. (3). Secondly, since these 
new oscillations of the magnetoresistance have been observed in weak magnetic fields and 
for weak modulations, we attempt a classical evaluation of the correction to the unperturbed 
energy En by the modulation Vx cos KxX + Vy cos Kyy using the classical equations of motion 
x{t) = xq + RcSm.{ijJct + ip),y(t) = yo + RcCos^uJ. + ip); Xq and yo are the classical center 
coordinates, Rc is the cyclotron radius, Uc = \e\B/m*, and (f is a phase factor. Without 
loss of generality we may take (f = 0. Then, if T is the period of the cyclotron motion, a 
straightforward evaluation gives 

fT/2 

<U> = (1/T) / dt[Vx cos Kxx{t) + Vy cos Kyy{t)] 

J-T/2 

= Vx Jo (KxRc) cos KxXq + Vy Jo (KyRc) cos Kyyo , (7) 

where Jo{x) is the Bessel function of order zero. In the weak magnetic field limit K^Rc ^ 1, 
Eq. (7) reduces to Eq. (3) for large n, i.e., for weak B. It is obvious that these three 
approximate ways of deriving the energy spectrum do not "see" its fine structure resulting 
from an exact numerical evaluation of the finite-difference Eq. (2) for a = i/j. Therefore, 
they are applicable if the corresponding small gaps are closed due to disorder. 



B. Hexagonal symmetry 

We assume that the Hamiltonian is given by 



ijo = {p + eAf /2m* + Vx cos KxX cos Kyy + Vy{l + cos2Kyy)/2. (8) 

.^OT Vx = Vy = Vo this reduces to the model studied experimentally by Fang and Stiles 
3]. If X and y are interchanged the energy spectrum, with Ky = 27r/a and Kx = 27r/v^a, 
of the corresponding tight-binding equation has been studied numerically, for all values of 



the magnetic field, by Claro and Wannier 11[ and has the same structure as that of the 
square symmetry. Here, in line with the case of rectangular symmetry, we assume that the 
small gaps of the energy spectrum are closed due to disorder and use again the tight-binding 
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description of Sec. II A. Corresponding to Eq. (2) we now obtain 

]^V^Fn{u^ + Mj^)[cos(27rpa + 7)^^+1 + cos(27rpQ; - 7)^p-i] + -^^n(4My)[^p+2 + 

^{E-E,,-\Vy)A,, (9) 

where 7 = Kx(?{ky + Ky/2). When a is integer Eq. (9) has the solution Ap = yloexp(i^p) 
with Aq and ^ given in A and the eigenvalue E is given by 

Enk^ = En+ ^Vy + V^Fn{u^ + Uy) COS 7 COS ^ + ^VyFn{4Uy) COS 2^, (10) 

where ^ — £^Kykx- We notice that for — 27r/a and Ky — 2^T^/3a, i.e. the usual hexagonal 
modulation, we have Fn{ux + Uy) — Fn{Auy) — F„(87r^£^/3a^). 

As in the rectangular case, we see that the Landau levels have broadened into bands with 
a bandwidth equal to {2Vx\Fn{ux + Uy) + Vy\Fn{'iUy)\) that oscillates with magnetic field and 
(large) n. Again the mean velocities Vx and Vy are finite 

Vx = -{VyeKy/h)Fn{Auy) siu 2^ " {Vxi'' K y / H) Fn{ux + Uy) COS 7 COS (11) 



% = -{Vxf Kx/h)Fn{ux + Uy) sin sin (12) 
this has important consequences for transport and will be detailed in the next section. 

C. The density of states 

The energy spectra given by Eqs. (3) and (10) are qualitatively different from the un- 
modulated spectrum, given by En, and from the corresponding ID modulation spectrum 
given by En + Fn{ux) cosK^xq. These differences are also reflected in the density of states 
(DOS) defined by D{E) — '^^nky^^iE — Enky^)- For a 2D modulation with rectangular 
symmetry, corresponding to Eq. (3), the DOS becomes 

oo /■27r 

D{E) ^DoJ^ d^{[VxFn{ux)f -[E-En- VyFn{uy) cos^f}-^/\ (13) 

n=0 

while for the one with hexagonal symmetry, corresponding to Eq. (10), the DOS is given by 



D{E)^DoJ2j^ dd[VxFn{ux + Uy)cos^f-i^E-En-f[l + Fn{Auy)cos2^]j 



-1/2 



(14) 



8 




FIG. 3: Density of states versus energy for Vx = Vy = 0.5 meV, ax = 800 A with ay = 800 A and 
ay = 1600 A for the sohd and dotted curves, respectively. The dash-dotted curve is the result for a 
ID modulation along the x-direction with the same period and modulation strength as in the 2D 
case. The magnetic field is B = 0.64 T. 

where Dq = LyLx/Tc^i"^. The quantities within the curly brackets in Eqs. (13) and (14) must 
be positive. 

In Fig. El we plot the DOS, given by Eq. (13), for various values of the parameters 
ttx, ay, Vx, and Vy. For comparison we also show the DOS (dash-dotted curve) corresponding 
to the ID modulation. The latter exhibits van Hove singularities at the edges of each Landau 
level (band) reflecting the ID nature of the electron motion in this band, since Vx ^ while 
Vy = 0, cf. Ref. 2 and 3. This is not the case for the 2D modulation: the electron motion is 
two-dimensional, since both Vx and Vy are different from zero, cf. Eqs. (5) and (6). That is, 
in the 2D case the DOS is finite, see also Ref. 4 b). As shown there, the DOS is qualitatively 
the same as the one shown in Fig. 01 if the periods are the same and the strengths are varied. 
This can be immediately deduced from the factors \4.F„(u^.) and VyFn{uy) that appear in 
Eqs. (3) and (11). The DOS for the hexagonal modulation is not shown since it's similar to 
the one shown in Fig. El 

In Fig. m we compare the DOS obtained from the exact energy spectrum with that 



9 




E (meV) E (meV) 



FIG. 4: (a) Density of states versus energy for a = 2/3 and energy level width F = (thin 
curves) and F = IK (thick curves). The solid and dotted curves are, respectively, the exact and 
approximate results. The DOS for a = 1 in (b) and a = 1/2 in (c) is plotted, respectively, for 
F = I.IK and F = 1.5K obtained from F = (eh/m*)\/B/TTii. The parameters used are the same 
as those in Fig. [21 

obtained using Eq. (3). In this comparison we include a level broadening by replacing 6{E) 
in the definition of D{E) by TfT/lE"^ + F^). In Fig. Efa) we show the influence of the level 
broadening on the DOS for a = 2/3 and different values of F specified in the caption. As can 
be seen, the subband structure disappears with increasing F and the exact and approximate 
result approach each other. That is, the gaps between the minibands in each Landau level 
are closed with increasing level broadening. Notice that this happens for quite small values of 
F compared to the cyclotron energy which is about 1 meV in this example. As shown in (c), 
the same behavior of the DOS occurs for a = 1/2. Notice also that this closeness between 
the exact and approximate DOS occurs despite the drastic difference in the corresponding 
energy spectra shown in Fig. 2. In addition, as shown in (b), for integer a the exact and 
approximate results for the DOS are identical. 
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III. TRANSPORT COEFFICIENTS 



A. Basic expressions 

For weak electric fields E, i.e. for linear responses, and weak scattering potentials the 
conductivity tensor a^i,{uj) in the one-electron approximation, has been evaluated in detail 
in Ref. 3: cr^^^uj) = crf^^iuj) + cr^^(c<j), /i, z/ = x,y. The contribution (J'^iy{uj) stems from 
the diagonal part of the density operator p. In a suitable basis < >= Tr^p'^J^j) = 
a'lj^^E^, where is the current density, and ct^^(co') comes from the nondiagonal part of 
p{p = p'^ + p'^'^). In general 0"^,^(co') = a'j^J {u) + cr^°'(u;), where a'j^J {oj) indicates diffusive 
contributions and cr™'(co') coUisional contributions. For the diffusive contribution we have 

'^^(0) = ^$^/c(i-/cMi5;cH^5, (15) 

c 



provided that the scattering is elastic or quasielastic, and for the collisional one 

<(0) = ^ E /c(l - k'W,,'{ai - a<:r, (16) 

for both elastic (Z,^ = /^/) and inelastic {f(^ ^ f(^i) scattering. W,^,^' is the transition rate 
between the unperturbed one-electron states | ( > and | C' >, ^ the volume of the system, e 
the electron charge, t{E() the relaxation time, and =< C,\r^ \ Q > the mean value of the 
/i-component of the position operator when the electron is in state | C > and has velocity 

=< C I "^At I C >• Equation (15) describes transport through extended states whereas Eq. 
(16) deals with transport through localized states and is absent in semiclassical treatments. 

The nondiagonal contribution o""^(co') to the conductivity is given by 



J2fd^-k)<(\^,\C><C\v.\(> 



1 _ (,P{Ec-E^,) 

X — 1=, 1=, lim-; T -. (17) 

Ec_ - E(^, e^o Ec_ - Ec_i + huj + le ^ ^ 

If we use the identity /^(l — fc^i) exp [[3{Ec^ — -E^)] = /c'(l ^ /c)' 1^^- (!''') takes the form of 
the well-known Kubo-Greenwood formula. 

Apart from their use in Ref. 3 for the ID modulation case, the above formulas have 
also been succesfully applied to various situations of electronic transport, such as hopping 
conduction [l^a, Aharonov-Bohm effect Qlb, quantum Hall effect |l^c, etc. 
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The resistivity tensor p^^y is given in terms of the conductivity tensor p = a ^ . We will 
use the standard expressions pxx — '-^yy/^iPyy ~ ^xx 

/S, and Pyx = -pxy = -(Jyx/S with 



^ ^xx^yy ^xy^yx' 



B. Analytical evaluations 

The scattering mechanism enters the conductivity expressions (15) and (16) through the 
relaxation time t{E() and the transition rate Wc^c, respectively; in contrast, Eq. (17) is 
independent of the scattering when the latter is weak 13[. 

We assume that the electrons are scattered elastically by randomly distributed im- 
purities. This type of scattering is dominant at the low temperatures of the reported 
experiments. Further, we expand the impurity potential in Fourier components, i.e., 
U{v — R) = f/q exp [iq.(r — R)], with ?7q = 27re^/e(g^ + klY^"^ corresponding to the 
screened impurity potential U{v) = (e^/er) exp(— fc^r); r and R are the electron and impu- 
rity positions, respectively, q = qxX + qyi), e is the dielectric constant, and kg the screening 
wave vector. 

Diffusive contribution. For weak modulation potentials Vx and Vy, which is pertinent to 
most of the reported experiments, we may use | ( >=| n,ky,C, > to evaluate the velocity 
matrix elements appearing in Eq. (15); the latter are given by Eqs. (5) and (6) for rectangular 
symmetry and by Eqs. (11) and (12) for hexagonal symmetry. As for the relaxation time 
t{E(^), it is defined by 1/t{E^) = W^CC'('^C ~ Though the Landau levels broaden 

into bands, this definition fails at the fiat-band conditions, when = f ^ = 0. For this reason 
we estimate it from the lifetime given by l/r(£'^) = In the limit kg ^ q, we obtain 

r = r(£'^) ^ {ni'^ff/NjUQy^'^ where Nj is the 2D impurity density and Uo ^ 2Tie^/eks. 
However, at weak magnetic fields we may use r as constant and estimate it from the zero- 
field mobility p: r = tq = pm*/e. 

We now use Eqs. (5), (6), and (15) with (Ly/n) Jq""^^^ dky and 

{Lx/nax) d^- The result for a^^J is 

^ J^x^''^'^'' ^ ^^"^""^^^^ r dkyUyd^ - U^^) sin^ e (18) 

The component aylf is given by Eq. (18) with x and y interchanged, and sin^^ replaced 
by sm^{i'^Kxky) = sin^(i^'a,a;o). In the limit of vanishing Vy Eq. (18) gives the result of 
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a ID modulation, afj = 0. If we neglect the weak ky- and ^-dependence of the factor 
fnkyi{^ — fnky^) wc obtain the simplified expression 

n 

The corresponding expression for (jyy^ is given by Eq. (19) with x and ?/ interchanged. 

For hexagonal symmetry we use Kx = 27r/a and Ky = 2tx j \f?)a. The results for a'^^ and 
a'^^ are similar to Eq. (18) and can be easily obtained using Eqs. (11) and (12 ) for the 
velocities. The result for a^^-^, corresponding to Eq. (19), is given by Eq. (19) with K^Vy 
replaced by K^V^ /2 and that for afj by 

where u = Sn^f/Sa'^. 

Collisional contribution. To evaluate this contribution to order V^^ we must use the 
perturbed wave function to order V^. The procedure for evaluating Eq. (16) is identical 
with that corresponding to the ID modulation detailed previously 3[. We have again < C, \ 
X I C > — < C I a; I C' >= ^^{ky — ky); the only new ingredient are the following matrix 
elements 

< nky^ I y I nky^ >= -^/Ky (21) 

and 



< nky^ I e^^- I n'k'y^ >\'= {n\/n'\)u- ""e"'^ ""(«) (22) 



where m = £^(g^ + g^)/2 and Cy = i'^Ky. 



We now use Eqs. (16), (20)-(21), and the standard expression for the transition rate 



q 

fj2 . 



We use the spectrum (3) and shift the argument of the cosines by £ K^Ky = 27ra, a integer, 
in the 6 function as well as in the factor fnky^i^ — fn'k'y^')- Then Eq. (16) takes the form 

l2 



^ h Aa^ 



poo 

J2 / rfMe-V'-"+^[L::'-"(M)]^ 

/ -'0 



/ dC, dky fnkyd^ " fn' ,ky+Ky+qy,(,-Cy(K^+qSi 

Jo Jo 



X 



X ^{Enky^ - En',ky+Ky+qy,i-Cy{K^+q^))- (24) 
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We proceed as follows. For weak magnetic fields involved in the problem the Landau- 
level index n is large and the major contributions to the sum over n' come from n' val- 
ues close n. With the asymptotic expansion of the Laguerre polynomials, e~"/^L„(M) a; 
(n^nu)~^^* cos{2^/ml — vr/4), it's an excellent approximation to take Fniu^j) ~ Fn'{u^). 
Then the 6 function becomes 

2Fn{ux)Vx sin Cx{Ky - qy/2) sin Cx{ky + Ky - qy/2) 
+ 2Fr,{uy)VySmCy{K^- q^/2)s\ucy{k^ + K^- q^/2)\. (25) 

The shift by i'^K^Ky — 27ra, a integer, in Eq. (25) and in the factor fnky^i^ — fn'k'^^') was 
made to stress the formal vahdity of Eqs. (24) and (25) for a integer. If we don't make 
it, we must put ^ Ky ^ in the sine factors and change En' ,ky+Ky+qy,e,-cy{K^+q^) to 
En',ky+qy,^-cyqj: wherever it appears. For a close to an integer though, one can reinstate 
and Ky in Eqs. (24) and (25) as shown. 

We now remark that the largest contribution to the integral over u in Eq. (24) comes 
from very small values of q^ and qy due to the factor exp(—u) or the factor l/V n'^nu in the 
asymptotic expression e~"[L„(ii)]^ ^ cos^(2-v/mZ — n / A) / Vnhiu. In addition, for the usual 
2D systems we have ks ~ 10^ /m which is much larger than these small values of qx and qy . 
With that in mind and in order to reduce the numerical work, we replace the S function (25) 
by a Lorentzian of width F and neglect in it and in the factor fnkyd'^-fn',ky+Ky+qy,^-cy(K^+ga:)) 
the terms (x q^ or oc qy. Alternatively, we may expand the S unction in powers of q^ and 
qy] then by far the leading contribution comes from the zero-order term given by Eq. (25) 
with qx — Qy — 0. In addition, we neglect the term q^ in U^. Further, from the sum over 
n' we consider only the terms n' — n and n' — n ± 1; the term n' — n gives the dominant 
contribution, about 90%. Then the integral over u can be evaluated and Eq. (24) takes the 
form 

X [Dn,n + (n + l)L>„,n+l + nL>n,n-l]} (26) 

where 

Dn,n' = fukyCi^ - fn',ky+Ky,^-CyKj/[{Enky£. - En' ,ky+Ky,e,-CyKj^ + F^] (27) 
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FIG. 5: Scattering between two cyclotron orbits which encircle one unit cell. 



As a result, when l-nP'jaxCiy = ^o/^ is an integer, the second and third terms in the 
argument of the 6 function in Eq. (25) vanish and entail n = n', i. e., the response is 
strongest when one flux quantum passes through an integral number of cells as observed 
In this case the factor [{...f + V^] in Eq. (27) becomes T^. 

A qualitative understanding of the enhancement of the coUisional conductivity for integer 
a = ^o/^ is as follows. In this case only scattering between the states \ky,^) and \ky + 

+ 1y^^ + cy{Kx, q^)) is allowed, cf. Eq. (24). These states correspond to cyclotron orbits 
separated by a distance aa^ which is a multiple of the lattice period. One example is shown 
in Fig. 5 for two orbits that encircle a unit cell. As shown, the orbits are in the same relative 
position with respect to the modulation lattice and correspond to electron states of the same 
energy. Since impurity scattering is an elastic process that leads to hopping between states 
of the same energy, the hopping between such cyclotron orbits for integer a contributes the 
most to the conduction and enhances the collisional conductivity. On the other hand, for 
a not an integer the position of the two orbits involved in the scattering process relative to 
the modulation lattice changes; accordingly the enhancement mentioned above is weakened. 

For those values of the magnetic field for which Fn{ux) vanishes we use the same wave 
functions and the spectrum (4). If the modulation periods are the same, we have Fn{ux) = 
Fn{uy) = 0, n — > n', and Eq. (26) holds with Dn,n±i 0. If the modulation periods are 
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not equal or if $o/^ is not an integer, Eqs. (25) and (26) hold only approximately. With 
all that in mind, the assumption that the small gaps are closed due to disorder, and for 
computational convenience, we use Eq. (26) as an approximation for all fields. 

For the hexagonal modulation we obtain again Eq. (26) but now the energy spectrum is 
given by Eq. (10). Further, and Uy are replaced by a and u = Sn^i"^ /Sa"^, respectively. 
For a^°J the result is given by Eq. (26) with ax replaced by a; u remains the same. 

The Hall conductivity. The evaluation of Eq. (17) for u; = is readily performed with 
the states | nk^^ > and the energy levels given by Eqs. (3) or (10). The only difference with 
the previous [2| calculation is that a factor exp[i{ky — k'y) / Ky]5^^^' appears on the rhs of Eq. 
(17) of Ref. 3 now written as < nky^ \ \ nk'^^' >. For rectangular symmetry we obtain 

^- = 17— + 1) / dky / fnky, U^iky, ^28) 

h Tia^ ^ Jo Jo [I + Xnx COS{KxXo) + Any cosQ^ 

where = V^e~'^''^'^L^l^-^^{u^) /fkOc, fJ' = x,y. We notice that for K/ = we obtain the 
previous ID result Q]. We also remark that Eq. (28) is valid for hexagonal symmetry with 
a^c ^ a, Ux = Uy = u = 87r^£^/3a^ and of course the different energy levels (Eq. (10)) that 
enter the factor fnky^- 



IV. NUMERICAL RESULTS 



We now present results for the various resistivity components p^i, using the standard 
expressions given at the end of Sec. III. A, and evaluating numerically the conductivities 
given by Eqs. (18), (26), and (27). For Figs. 3-7, and 10 we use the parameters of Ref. 8. 
They are: electron density Us = 4.5 x 10^^/m^, temperature T = 5.5 K or T = 1.6 K, = 
ay = 804A, and mobility /i = 70 m^/Vs. The corresponding parameters for Figs. 8-9, taken 
from Ref. 5, are = 5.1 x 10^^/m^, temperature T = 4.2 K, ax = ay = 2820A, and mobility 
/i = 140 m^/Vs. The relaxation time at zero magnetic field tq is estimated from the sample 
mobility as Tq = m */i/e. Then the level width is F = {eBNiU^ /Tihy/^ = {eh/m*)^ B /up. 

In Figs. 3 we plot pxx and pyy as function of the magnetic field B with Vx = 0.5 meV for 
constant r = tq. As indicated, the various curves correspond to different Vy and the dotted 
one, marked ID, represents the ID limit obtained with Vy = 0. The prominent peaks in the 
2D case, marked by the integral value of a, result from the coUisional contribution to the con- 
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B(T) 



FIG. 6: Resistivity components pxx sjnd Pyy B-S cl function of the magnetic field B for fixed Vx = 0-5 
meV and varying Vy. The dotted curve is the ID hmit obtained with Vy = 0. The prominent peaks 
in the 2D case are marked by the integral value of a = ^o/^- 

ductivity. The smaller peaks, between these values of a, correspond to the commensurability 
or Weiss oscillations. Notice how the prominent peaks of pxx, in the 2D case, remain rather 
insensitive to changes in Vy-. this is so because they result from the collisional contribution 
^coi Tf^\^[Q\^ depends very weakly on Vy through the energy spectrum. The apparently drastic 
difference between the two figures results from the fact that a^^ ^ axy makes S change little 
and Pxx = o'yy/S while Pyy = cTxx/S. Upon reducing Vy the contribution afj ~ Vy, given 
by Eq. (18), is affected drastically whereas 0"^°' is not. 

In Figs. 4 and 5 we plot again pxx and Pyy as function of the field B with Vx = Vy = 0.5 
meV, for the same r and F oc B^^"^ as in Fig. [Ul but with the period ay being doubled 
from panel to panel as indicated. The solid curves give the total resistivity, the dashed 
ones the diffusive contribution, defined by p^*f = a'^l^ /S, and the dotted ones the collisional 
contribution, defined by p^'^' = a'^^ / S. Notice how the prominent peaks move to lower fields 
with increasing ay as explained after Eq. (26); in panel (d) they have disappeared. The ID 
limit shown in panel (d) is obtained with Vy = Q and the difference in the B dependence 
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FIG. 7: Resistivity component 

Pxx ^ 

function of the magnetic field B for fixed l/j.. = = 0.5 
meV and different periods ay as indicated. The sohd curves give the total resistivity, the dotted 
ones the diffussive contribution, and the dash-dotted ones the collisional contribution. The thinner 
curves in panel (d) are for the ID limit {Vy = 0.) The prominent peaks in the 2D case are marked 
by the integral values of a. 
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FIG. 8: The same as in Fig. [3 but for the resistivity component pyy. 



between p^^ and pyy is related to that of the corresponding conductivity contributions. One 
of them, afj given Eq. (18), is affected drastically by changing Vy and/or the period ay 
which enters the factor sin^ ^, the others very weakly. 

We now look more closely at the experimental results of Refs. 8 and 5 with the parameter 
sets specified above. The parameters Vx and Vy are not known. In Fig. IHlwe have shown 
the total resistivity for \4 = K = 0.5 meV, constant r = tq. In Fig. El we plot again 



versus B but now we show the contributions pfj and 



p™' as well. In addition, we 
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take Vx = Vy = 1 meV, and 1/r ~ F ~ B^/"^. Because af^J <^ a™', the difference in the 
total resistivity is very small between the two sets of modulation strengths. However, the 
oscillation amplitudes^ in are higher in the present case and p^x increases more slowly 
with B as observed 8]. Upon closer inspection we see that the prominent peaks, marked 
by the integral values of a, result entirely from the collisional contribution cr^°'. As can be 
seen in Fig. IHlof the next article Q|, the amplitudes and positions of these peaks agree well 
to very well with the experimental results. Notice also how the Weiss oscillations of p^^^ 
and between these peaks, are in antiphase. These experimental results for pyy and pxx 
are for two orthogonal crystal directions [Oil] and [Oil] taken from different samples. They 
have similar structures and the curves may be fitted theoretically using slightly different 
potentials. A direct comparison between experimental and theoretical results is made in 
Fig. ini of the next article jiM] for t = tq but qualitatively the agreement is the same for 
r oc Notice, however, that the theoretical oscillation amplitudes and the overall value 

of Pxx in the low-B region, below a = 2, agree less well with the experimental ones than 
those in the high-i? region. 

In Ref. results are given for temperature 1.6 K and otherwise the same parameters. 
We show the calculated pxx for this case in Fig. ^1 As can be seen, lowering the temperature 
makes visible all prominent peaks marked by arrows for a = 1, ...,8. Their positions occur 
at fields B = 0.64, 0.32, 0.21, 0.16, 0.13, 0.11, 0.09J).08 T and compare very well with the 
experimental ones, see Fig. IHlin the next article [lj|. To see more clearly the oscillations we 
replot, in Fig. ^2 the resistivities in the low-field region of Fig. ^Jas a function of 1/B. 

The temperature dependence of the oscillations is shown in Fig. [T21 The solid, dotted, 
and thin solid curves correspond to T = 5, 10, 20 K, respectively. As can be seen, these new 
oscillations are more robust than the Weiss oscillations and persist at T = 20 K. However, 
their damping with T is weaker than the observed one Q|. 

In Fig. ^Jwe plot pxx in the manner of Fig. El but for the parameters of Ref. 5 involving 
the much longer periods ax = dy = 2820 A. The modulation strengths are Vx = Vy = 0.2 
meV and very close to those used in Refs. 5 and 6. The agreement with the experimental 
2D results of Ref. 5 is very good: below approximately B = 0.5 T we have the Weiss 
oscillations and above it the Shubnikov-de Haas ones. One noticeable feature here is the 
absence of the prominent peaks for integral values of a. This is so because the much longer 
periods involved make a = 2TTi'^/axay integer for much smaller values of B. For instance. 
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B(T) 

FIG. 9: Resistivity component 

Pxx cl function of the magnetic field B for fixed Vx = Vy = 1 
meV and T = 5.5 K. The sohd curves give the total resistivity, the dotted ones the diffussive 
contribution, and the dash-dotted ones the collisional contribution. The prominent peaks are 
marked by the integral values of a = <I>o/<5. 
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FIG. 11: Resistivity component 

Pxx ^ 

function of inverse magnetic field 1/B. The curves are 
marked as in Fig. IIUI The inset shows the peak position versus \/B. 
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FIG. 12: Resistivity component 

Pxx ^ 

function of the magnetic field B. The solid, dotted, and 
thin solid curves correspond to T = 5, 10, 20 K, respectively. 



21 



20 




I . ..>.-lw-^ ., .--T~, .-'-' l , l "'~T-,. l 

0.0 0.2 0.4 0.6 0.8 1.0 
B(T) 



FIG. 13: The same as in Fig. I^with the parameters of Ref. 5 and Vx = Vy = 0.2 meV. 

a = 1 occurs at B = 0.05 T and the corresponding peak is not resolved. The agreement is 
also as good if we use ax = ay = 3650 A and otherwise the same parameters pertaining to 
another sample. 

Reference 5 reported results also for ID modulations. As mentioned earlier, we can obtain 
the ID limit from the present 2D results by considering a vanishing Vy. In Fig. ^^we show 
the ID limit of the total pxx and pyy for a^,. = 2820 A, Vx = 0.5 meV, and Vy = 0. Although 
the agreement between theory and experiment is very good, it must be noticed that it is 
obtained with Vx = 0.5 meV and not Vx = 0.2 meV that we used in Fig. El Since the 
the ID or 2D modulations are produced by illumination of the samples, we expect them 
to have the same strength. If we use Vx = 0.2 meV we can obtain good agreement if we 
use a T smaller by about a factor of 2 in Eq. (18). As stated in Refs. 5 and 6, this may 
be an indication that in this very high mobility samples the fine structure of the energy 
spectrum, that the present theory neglects, is partially resolved. However, the experimental 
data was taken at T = 4.2 K and, as no such fine structure has been observed above mK 

Finally, in Fig. ^1 we show the Hall resistivity pyx for the parameters of Fig. IHl As in 
the case of ID modulations, it exhibits very weak oscillations. They are better seen in the 
inset which shows the derivative dpy^/dB versus B. The triangles on the x axis mark the 
positions of the integral values of a for which enhanced oscillations are observed. 
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FIG. 14: The ID resistivity component p^^ as a function of the magnetic field B obtained with 
ax = 2820 A, K = 0.5 meV, and Vy = 0. 
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FIG. 15: The Hall resistivity pyx versus magnetic field B with the parameters of Fig. IHl The inset 
shows the derivative dpyx/dB versus B. 

V. CONCLUDING REMARKS 

We presented a theory of magnetotransport in 2D superlattices using the energy spectrum 
and wave functions that result from the tight-binding difference equation when the parameter 
a = $0/"^ is an integer. As emphasized in the text and supported with the results for the 
DOS shown in Fig. 4, the description holds approximately for all fields if we assume that 
the small gaps in the energy spectrum are closed due to disorder. The reasonable-to-good 
agreement with the experimental results strongly supports this assumption. 

As detailed in the text, the prominent peaks, for a = $o/^ integer, result from the 
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collisional contribution to the conductivity 0"^° , require sufficiently s/iori periods, and depend 
very weakly on the value of the modulation strengths Vx and Vy. Upon increasing the period 
along one direction we showed how they move to lower fields. Accordingly, for periods 
between 3000 A and 4000 A these peaks occur at much smaller magnetic fields and are not 
resolved |5jj]. The agreement between our results and the experimental ones, as presented 
in Ref. 8 and detailed in the next article is good for the peak positions at all fields. 
The oscillation amplitudes agree well at relatively high fields but less well at low fields. As 
shown in Fig. [T21 these oscillations are quite robust with respect to the temperature but 
their damping with temperature is weaker than the observed one. 

Between the oscillations for a = $0/"^ integer we have the Weiss oscillations. The 
relative phase between those of p^x and those of pyy depends on the values of the modulation 
strengths, cf. Fig. IHlin which the period is the same for all curves, and of the modulation 
periods, cf. Figs. 4 and 5 in which the modulation strengths are the same for all panels. 
We notice that as Vy becomes smaller and smaller than K., the oscillations resemble more 

n 

closely those corresponding to ID weak modulations [1^, cf. Fig. IHl The results for the 
latter can be extracted from the present 2D ones if we take the modulation strength along 
one direction to be zero. 

The relative phase between pxx and pyy for ID and 2D modulations depends strongly on 
the ratio of the modulation strengths Vx and Vy. Since one resistivity component vanishes for 
ID modulations, this conclusion could not be reached by studying only the latter. Similar 
results were reported in Ref. 6. 
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